Exact solution of a linear molecular motor model driven by two-step fluctuations and 

subject to protein friction 
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We investigate by analytical means the stochastic equations of motion of a linear molecular motor 
model based on the concept of protein friction. Solving the coupled Langevin equations originally 
proposed by Mogilner et al. (A. Mogilner et al., Phys. Lett. 237, 297 (1998)), and averaging over 
both the two-step internal conformational fluctuations and the thermal noise, we present explicit, 
ryj ' analytical expressions for the average motion and the velocity-force relationship. Our results allow 

for a direct interpretation of details of this motor model which are not readily accessible from 
numerical solutions. In particular, we find that the model is able to predict physiologically reasonable 
values for the load-free motor velocity and the motor mobility. 

', PACS numbers: 87.16.Nn, 02.50.-r, 05.40.-a, 87.10. +c 
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I. INTRODUCTION 

> 

There is currently widespread interest in molecular motors, both from a biochemical-physiological and a physics 
point of view. Whereas the former is mostly concerned with the molecular structure of motors and their structural 
interplay with the support on which they move, physicists study the non-equilibrium transport properties of motors 
and their physical interactions with the support, such as load- velocity relations or adhesion forces between motor 
• and support. Molecular motors, in general, are energy-consuming, non-equilibrium nanoscale engines, which are 
encountered in various dynamical processes on the intra- and intercellular level 0, Q ; for a recent review of the more 
physical aspects see, for instance, Ref. 0- Such motors are responsible for intracellular transport of molecules and 
small vesicles in eukaryotic cells; they are powering genomic transcription and translation, cell division (mitosis), or 
the packaging of viral DNA into nanoscale transport containers (capsids) 0, HI 0, Gi @j El ■ Larger assemblies of motors 
working in unison are responsible for the motility of, e.g., bacteria, the y p lay a role in cell growth, and they are 
responsible for muscle contraction leading to macroscopic motion 

Linear motor proteins like myosin, kinesin, dynein, or DNA helicase or RNA polymerase are driven by the cyclical 
hydrolysis of ATP into ADP and inorganic phosphate and wander along linear, polar biomolecular tracks such as 
actin filaments, microtubules, RNA or DNA. The motion is typically associated with two or multi-step conformational 
changes in the motor protein in interaction with ATP/ ADP and the filament /support, and takes place in a thermal 
environment subject to viscous forces. 

Modern experimental techniques in biology and biophysics, in particular single biomolecule manipulation by for 
example optical tweezers or micro-needles, and single particle tracking methods , ha ve yielded considerable insight 
into the mechanism and the relevant physical scales in molecular motor behavior [llL Il2l Iba Hfl Ha lL7l fl8j| . The 
typical size of a molecular motor is of order 10 - 20 nm, moving with a step size of order 8 nm, e.g., kinesin on 
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microtubules, with one ATP molecule hydrolyzed on the average per step. The velocities of molecular motors range 
from nm/sec to /im/s and the maximum load is of the order of several pN (e.g., ~6pN for kinesin on microtubules). 
However, the latter can reach up to 57pN for the rotating packaging motor of bacteriophages [ljj. The time scale of 
the chemical cycle is a few ms and the average energy input from the ATP-ADP cycle of order 15-20 kT. 

A biomolecular motor represents an interesting and ubiquitous non-equilibrium system operating in the classical 
regime and is thus directly amenable to an analysis using standard methods within non-equilibrium statistical physics. 
Physical modelling of molecular motors has thus been studied intensively in recent years both from the point of 
view of the fundamental underlying physical principles and with regard to specific modelling of concrete motors 
[8l I20I l2ll 12^ . 12^ . [2^ . l25l 12^. l27 l [28 l l29j . More recently, the concerted action of multiple motors has been considered, 
such as the action of elastically |31j and rigidly 01 coupled motors, for instance, in muscles . Motors interacting 
with freely polymerizing microtubules or actin filaments give rise to rich pattern formation such as asters |33|. and 
are responsible for the formation of the contractile ring emerging during cell division 0, [3j| . 

The most common statistical approach to molecular motors is that of a ratchet model 0, |24| |. mimicking the 
periodically alternating energy landscape (given by the interaction potential with its support) perceived by the motor 
during its mechanochemical cycle. Such ratchet models date back to Smoluchowski |3fi| and Feynman [37j |. and 
Huxley's pioneering work |38j | on motor proteins actually corresponds to a Brownian ratchet We note that 
ratchets play a much more general role, and real-space ratchets may even be used on the microscale for particle 
separation |39ll40| . 

An alternative motor model can be based on protein friction 0, 0, Ell E2I l43l |44| . This concept relies on the 
idea that due to the weak chemical bonds forming between motor protein and the polar actin or microtubule track, 
after elimination of the detailed degrees of freedom, an effective friction ( p builds up between motor and track. This 
protein friction £ p acts like a linear friction if the associated time scale of motion is longer than the characteristic 
time of the kinetics of motor-track bonds. If not, no protein friction can build up, and the motor is only subject 
to the smaller viscous drag due to the environment. The protein friction is therefore highly non-linear. On the 
basis of this scenario, Mogilner et al. recently studied a simple two-step linear molecular motor represented 

by two coupled overdamped oscillators driven by a two-step Markov process alternating between a relaxed and a 
strained state of the oscillators and embedded in a thermal environment represented by additive white noise. The two 
subprocesses are associated with internal conformational changes of the motor protein such that one subprocess is 
slow, allowing protein friction to be established, while the other is fast and only subject to solvent friction. By means 
of a numerical analysis Mogilner et al. show that the system acts like a motor and can carry a load. However, unlike 
the ratchet models, which operate with an attachment to a periodic polar protein filament, the model of Mogilner et 
al. only needs a 'passive' groove in order to perform directed motion, and the 'ratcheting' comes about by assuming 
the asymmetric internal velocity fluctuations, which are then rectified by protein friction. In that sense, it is a robotic 
model of molecular motors. 

In the present paper we reanalyze the motor model by Mogilner et al. from a purely analytical point of view 
and derive explicit expressions for the motion of the motor and the velocity-load relationship. Using the biological 
parameter values quoted by Mogilner et al. we show that the model gives rise to physiologically reasonable values for 
the motor velocity, whereas our analysis leads to a correction of the maximum load force by an order of magnitude 
m comparison with the numerical results obtained in Ref. 0. This discrepancy is associated with a difference in 
the dynamics of the analytical model as compared with the numerical simulation. Allowing for a larger relaxation 
rate the analytical result for the maximum load force approaches the biological regime. The paper is organized in 
the following manner. In Sec. [H] we introduce the model. In Sec. IIIII we solve the model analytically. In Sec. II VI we 
discuss the results and compare with Ref. |45j| . The paper ends with a summary and a conclusion in Sec. IVl 



II. MODEL 

The motor model based on protein friction which was introduced in Ref. |45|. is defined as follows, compare Fig.^ 
Assume that the mechanochemical cycle of the motor protein walking along the track made up of an actin filament 
or microtubule can be pinned down to the periodical switching between two states, and that each of these two states 
can be described by two motor heads connected by an effective spring representing the backbone of the motor protein 
between these heads. From the strained state (S), characterized by a rest length L s , the motor protein converges 
towards a relaxed state (R) with rest length L r > L s , i.e., the distance between the motor heads increases. This 
process is slow enough to make sure that the adhesion between the motor's 'working head' (white circle in Fig. ^) 
and track stays intact, in such a manner that asymmetric motion with respect to the track is achieved (stick). In 
contrast, during the fast 'power stroke' from the relaxed state back to the strained state after hydrolysis of ATP, the 
protein friction is broken and both heads move under the low Reynolds number conditions of the environment, such 
that both heads symmetrically approach each other and assume the rest length L s (slip) . 
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This model can be cast into the two coupled Langevin equations: 

C(t)i(t) = -f + k(t) y{t) ~ X ^~ m +N x (t), (2.1) 

CvVW = -«(*) 2 ^i/Wi ( 2 - 2 ) 

in which we have introduced the time-dependent friction coefficient £(t) in comparison to Ref. |45j for convenience, 
to account for the cyclical attachment to the track. In Eqs. (|2.1I2.2II the variable x and y represent the positions 
of the two heads of the motor molecule along the track, corresponding to the equations of motion of two coupled, 
overdamped oscillators. The coordinate of the idle head y is associated with a viscous friction drag coefficient ( v 
of order 6irr)r (Stokes law), where rj is the viscosity of water and r the size of the motor protein head. The same 
friction acts on the working head during the fast conformational change R— >S, whereas during the slow process S— >R, 
it experiences the protein friction drag with coefficient ( p [2ll l24l l4ll l42l l43l , corresponding to a stick-slip motion 
of the working head. The model equations (|2.1I2.2|) are driven by thermal noises N x (t) and N v (t), with {N X:V (t)) = 0, 
representing the ambient environment with correlations 

{N x ,y(t)N Xiy (t')) = 2k B TC P , v S(t - t% (2.3) 

balancing the friction terms by means of the fluctuation-dissipation theorem. We note that during the detached 
straining step, the role of working and idle heads may be exchanged (i.e. , the motor heads may be turned around a 
common axis), as was recently demonstrated for kinesin motor heads |30|. 

The conformational changes of the motor driven by the ATP-ADP hydrolytic cycle and the cyclical attachment 
to the substrate correspond to a continuous two-state Markov process for the time dependent rest length L(t), the 
time dependent spring constant k(t), and the time dependent protein friction C(i), alternating between state R with 
rest length L ri spring constant k r , and protein friction C = Cp an d state S with rest length L s , spring constant 
k s , and viscous friction ( = £„. The power stroke conformational transition R— >S driven by the ATP hydrolysis is 
characterized by the rate g s ; the relaxational conformation change S— >R has the rate g r . 

The relevant biological parameters quoted in Ref. J45j, entering Eqs. i|2-l[) and (|2.2() . are: Rate of hydrolysis 
g s ~ 10 3 s _1 , relaxation rate g r ~ 10 3 s _1 , spring coefficient in relaxed state k r ~ 0.01 pN/nm, spring coefficient in 
strained state k s ~ 0.5 pN/nm, rest length in relaxed state L r ~ 40 nm, rest length in strained state L s ~ 20 nm, 
viscous drag coefficient ( v ~ 10~ 6 pN s/nm, protein friction drag coefficient ( p ~ 5 x 10~ 5 pN s/nm, and load force 
f ~ ±1 pN. For further discussion of the model and parameter choices under biological conditions we refer to ref. 
|45| . An important difference between the model of Mogilner at al. and the present one in Eqs. (|2.1|) and l|2.2|l is that 
we take £(t) = ( v during the entire duration of the strained state S, while Mogilner assumes = ( v only in the 
first short time interval as the spring contracts (a time interval of order t s n p ~ Q v /k s \ in the simulations of Ref. |45| 
t s iip is taken infinitcsimally small), after which the protein bonds will form and protein friction take over, = (p. 
The two models will be similar if the relaxation time g^ 1 is of the order of t s ii p . 

III. ANALYSIS 

In this Section, we present a solution scheme for this motor model. The results obtained arc then further analyzed 
in the following Section. 

A. Analytical solution 

The motor equations l|2.1() and 12.211 are readily analyzed by (i) solving Eq. I|2.1|l for y(t) and deriving dy/dt, 
(ii) eliminating y in Eq. (|2.2|) and setting the two expressions equal to one another. We thus obtain the following 
equations for v x = dx/dt and v y = dy/dt: 

HCCv/k)v x + (C + Cv~ 2(Uk/k 2 ) + 2CCv/k)v x = 

-/(l - 2Q v {k/k 2 )) - ( V L + 2( v N x /k + (1 - 2Q v k/k 2 )N x + N y , (3.1) 
CvV y = -f-(v x +N x + N y . (3.2) 

Denoting the initial velocity at time t — by v x q Eq. I|3.1|l is readily solved by quadrature and together with Eq. 
(|3.2() we obtain 



v x (t) = 



k(t) 



-lit) 



Vx0 W) + CvL{t ' ] ^ ( *' )]e7( ° 



(3.3) 
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v y (t) = ~~ + h N *® +N vW ( 3 - 4 ) 

which form the basis for our discussion. 

We have introduced the renormalized load force /, the renormalized noise N, and the integrated spring and friction 
constant 7: 

/=/(l-2C„(fe/* 2 )), (3.5) 

N = 2( v N x /k + (1 - 2( v k/k 2 )N x + N y , (3.6) 

7 (t) = (1/2&,) / fc W + / k{t')/2{{t')dt'. (3.7) 
Jo Jo 



B. General properties 



We note various general features of this solution. First, both the load force / and the thermal noises N x and N y are 
renormalized by the fluctuating spring constant k. Secondly, the thermal noise basically enters additively and entails 
thermal fluctuations of the velocities. Since the stochastic conformational changes giving rise to the fluctuations of 
k, L, and £ are independent of the thermal fluctuations we can in order to monitor the time dependence of the mean 
motor velocity with impunity average over the noise. Note that the heat bath, of course, still enters through the 
friction coefficients. In the long time steady state limit we can ignore the initial terms and we obtain the reduced 
equations for the thermally averaged velocities 

«*(*) = ~ H 2^) Jo + C ^ ( * /)]e7( °' (3 ' 8) 

«»(*) = + fc(t y (t) f *'[/(*') + U(t')]e^ - f, (3.9) 

Z( ,v JO <*v 

which we proceed to discuss. 



C. Constant spring constant and rest length 



Let us first as an illustration consider the case of a constant spring length L, a constant spring constant k, and a 
constant protein friction ((t) = ( p . In this simple case 7 (i) = k((( p + ( v )/2( p ( v ))t and the load force is unrenormalized. 
We obtain 



/ 



Cp + Cv 



(3.10) 



Here the load / after a transient period drives the idle and working heads with a constant mean velocity. Defining 
the mobility according to 

«x = -A*/, (3.11) 
we infer the mobility in the absence of conformational fluctuations 

In the absence of a load for / = 0, the mean velocity vanishes and the system does not move, i.e., we do not have 
motor properties. This is also a statement of the second law of thermodynamics expressing the fact that we cannot 
extract work from a system in thermal equilibrium. In the case of constant fc, constant L, and constant Q the coupled 
Langevin equations describe the temporal fluctuations of a system in thermal equilibrium. The motor property is thus 
necessarily due to the fluctuating spring constant k(t) and fluctuating length L(t) characterizing the conformational 
fluctuations in combination with the cyclical attachment described by the fluctuating friction coefficient ((t). 
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D. Fluctuating spring constant, rest length, and protein friction 

The idea behind the model is that fluctuations of the spring constant k(t) = k s ,k r (k r < k s ) and rest length 
Lit) — L s ,L r [L r > L s ), modelling the ATD-ADP driven conformational changes, provide an energy source. In 
combination with the synchronized stick-slip mechanism modelled by a fluctuating protein friction £(t) = £ p , £ v 
(Cp > Cv), this process can drive the system in the absence of a force. This mechanism is modelled by the two-step 
Markovian process S<->R with relaxation rates g r for S— >R and g s for R— >S. The master equations for this process 
denoting the corresponding probabilities by P s {t) and P r (t) thus take the form: 

d ^M=g sPr (t)-g r P s (t) 1 (3.13) 

dP r {t) 
dt 

with stationary solutions 



g r P s {t)-g s Pr{t), (3.14) 



Ps = — J— , (3.15) 
9 s + g r 

P r = _£l_. (3.16) 

9s + g r 

The stationary mean value of, e.g., the spring constant, is thus given by 

fe st = k s P s + k r P r = fcsgs + fcrgr . (3.17) 

5s + 9r 

For the present purposes it turns out to be more convenient in discussing the conformational transitions to focus on 
the probability distributions P s (t) and P r (t) characterizing the residence of the system in either the strained state or 
the relaxed state at a time t. The distribution is exponential in time and we obtain properly normalized 



P s (t)=g r e- grt , (3.18) 
IUh /M :l '. (3.19) 



The mean values of the residence times are then given by 



9s 

The mean value of k may then be obtained as a time average: 



(*>. = -. (3-20) 

9r 

(t)r = —■ (3.21) 



v _ fcs (t) s + k r (t) r _ k s /g r + k r /g s 

W ~ (t)s + (t)r l/9r + l/ 9s ' ( ' 



in accordance with Eq. H3.17fl . 



E. The motor property without a load 

Here we establish the fundamental motor property of the model in the absence of a load. For / = we have from 
Eqs. and O, 

"° = -^%f i'^ 4 '' 67 '''' (3 ' 23) 

„ = + *(i)e-*" f'^nyf) (3M) 

z( ,v Jo 

At a superficial glance it looks like the motor heads move in opposite direction. However, subtracting the velocities 
and noting that 

Ht) k(t) 
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we obtain 



u° x -v° v = -y(t)e-^ / dt'L(t')e^ 



(3.26) 



Finally, assuming ergodicity (to be established later) and time averaging in combination with partial integrations we 
have for T — > oo 



(«2> - K> 



i 

rj 



i 



dtj^e-^ I dt'L(t')e^ t "> 
dt'L{t')e^ 



/ dt— 


e -7(*)" 


/o dt 





-- / L(t)dt = 0, 



(3.27) 



where the last step corresponds to an integration by parts (note that j(t) is monotonically increasing, and therefore 
r 1 ^ 

T _1 J Q L(t') exp{7(i') — 7(t)} < T -1 [£(t)] — > 0), and we conclude that the average velocities of the two heads 



are in fact identical: The working head and the idle head move together. 

Next we derive an explicit expression for (v®). Introducing the auxiliary fluctuating variable 



1 1 



and using the above result we obtain by adding Eqs. I|3.23ll and (|3.24|) : 



dt'{e-^-^a{t)L{t')). 



2 Jo 



(3.28) 



(3.29) 



Here (• • •) denotes an ensemble average with respect to the conformational fluctuations. Since L(t), k(t), and are 
governed by the same stochastic process, with the values L s , k s , and L r , £ p , k r in the strained and relaxed states, 
respectively, we obtain from the expressions for a and 7 in Eqs. i|3.28[) and (|3.7() 



(3.30) 
(3.31) 
(3.32) 







' 1 


1 " 


a r 


2 






a s 


= o, 






lr 


= T 


' 1 


1 " 


is 


k s 







(3.33) 



Denoting the jump times for the transitions between the strained and relaxed state by t n , n = 1, 2, • • • and assuming 
that the system is in a relaxed state at < t < t\ we have 



L(t) = (L r -L s )J2(-l) n 5(t-t n ), 



71=1 



and inserting in Eq. 1)3.29(1 we obtain 



<«2) 



f2(-lT(a(t)exp(-J%(t')dt^y 



(3.34) 



(3.35) 



where t/v < t < tjv+i- Note that the exponential term is just a more complicated way to write exp{— j(t) + 7(t n )}, 
which will be useful below. For even N the system is in the relaxed state for tjv < t' < t with probability P r . 
Similarly for odd N the motor ends in the strained state which occurs with probability P s . Introducing the time 
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interval r„ = t n+ \ — t n , noting that the residence distributions are statistically independent, and introducing the 
notation 



R = (exp (-7 r T)) r 



drP,.(r) exp (->r) = — 

/0 5s + 7r 

oo 

5 = (exp (-7„r))g = / rfrP s (r) exp (~7 s t) = 
Jo 



.9r + 7s ' 



the mean velocity can be expressed in terms of geometrical series, 

L r — L s 



(v° x ) = 



P r a r R{l - S) J2( SR ) n - p sa s S(l - R)J2( SR Y 



n=0 



n=0 



or summing the series (completing N — > oo), 

L r - L s P r a r R(l -S)- P s a s S(l - R) 
{Vx) 2 1-RS 

Inserting P r , P s , R, and S from Eqs. (|3~ll)|) . ESEJl, and (|3~37|) we arrive at 

2 (g s j s +g r j r + j r j s )(g s + g r )' 



(3.36) 
(3.37) 

(3.38) 

(3.39) 
(3.40) 



First we note that the expression vanishes for a(t) = j(t) thus corroborating the validity of the time average in Eq. 
H3.27(l and establishing ergodicity. Finally, inserting a r , a s , j r , and j s from Eqs. I|3.30(l . (|3.31(l . I(3.32|l . and l(3.33|l we 
obtain for the explicit expression for the motor velocity in the absence of a load 



(«2> 



(L r — L s )((^p — Cv)k r k s g r g s 



2(ffr + g s )[2g s ks(p(v + K{Cp + (v){g r (v + k„)] 



(3.41) 



F. The motor property with load 

We next turn to the case of a load force applied to the motor. First we establish that in the presence of the load 
the two heads of the motor move together with the same average velocity. From Eqs. I(3.8jl and 1(3. 9|1 we obtain 



v x v y = v° x - v° y - j^f^L f rff /(i')er(O + L. 

S>1> Jo S>i> 

Inserting / from Eq. i|3.5(l and averaging over time we have, using Eq. I|3.27(l : 

\v 1 JO Jo VM 



(3.42) 



(3.43) 



Performing the integrals by partial integration along the same lines as in the load-free case, the first two terms in Eq. 
(13.43(1 cancel and we find in the limit T — > oo, 



(V X ) = (Vy), 



(3.44) 



i.e., the two motor heads move together with the same average velocity. 

We now turn to the evaluation of the load-velocity relationship. From Eqs. (|3.8|l and (|3.9|) and inserting / we have 



Jr- I ill 



2C 



v Jo 



,-(7(*)-7(*')) 



..0 / , / 



v « = <-Z + xJo dt c 



,*(*) 



dt' \k(t') 

d (I 



-(7(*)-7(*')) 



(3.45) 
(3.46) 
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From the synchronization of the stochastic processes we obtain the identity 

11 l/k s -l/k. 



k(t) k r 7s - 7r 

and therefore 



(i(t) - 7r ), (3.47) 



dt k{t) 7 S - j r 

Inserting into Eqs. I|3.45[) and (|3.46|) and averaging 



d 1 l/k s — l/fc r .. . . .„ 

- ' i(t). (3.48) 



( v \ - ( v 0\ L(*® [ d t ' e -(7(*)-7(t')) V 

- f 1/fcs ~ 1/fcr /§4 T ^' e - (7(t) - 7(t,)) 7(^)^ , (3-49) 

7s - 7r \ CW Jo / 



-(7(*)-7(*')) 



K) = K) - T + wil 1 ! dt ' e 

+ f v k *-VK /m r* ^ e -(7w-7(*'))^( t ')\ 

7s - 7r \ Cu Jo / 
The first integral in Eqs. (|3.49(l and H3.50j) has the form 

/i = dte-^-^'A = ^(*)^ *exp ^- jf dt" 7 (t") 



(3.50) 



(3.51) 



and is performed by breaking up the integration over 7 in the exponents and averaging over the time segments yielding 
again a geometrical series in terms of SR. We obtain as an intermediate result 



^ 7r 7s / n=0 \ 7s 7r / n=Q 



(3.52) 



and performing the sum and inserting i? and S* from Eqs. (|3.36() and i|3.37[) . 



Zi = b r P r (g r +g s + 7s) + fr s P s (gr + ffs + ->) ^ ^ 

3s7s + ffr7r + 7r7s 



The second integral has the structure 

h = (c(t) J^dt'e-^-^^{t')\ , 
and was performed in the load-free case in Eqs. H3.35fl to (|3.40|) . We found 

T Us - j^PrPsigs + ffr)(Cr7s ~ C s j r ) 

H — : : ; — ; ■ 

9sls + grlr + Iris 



(3.54) 



(3.55) 



It is again convenient to introduce the mobility /j. according to the relation 

(v x ) = (v° x ) - M /, (3.56) 
and we obtain inserting c(t) = b(t) = k(t)/( p (t) for {v x ), or c(t) = b(t) = k(t)/( v for (v y ) 

{k r /2( v ( p )P r (g r + g s + 7,) + {k s /2$)P s (g r + 9a + 7r) 



/' : 



5s7s + 9rlr + jrjs 
P r Ps(g S + 9r){{K/Cp)js ~ (fcs/Ci»)7r)(l/fcs ~ l/fcr) 
5s7s + g r i r + i r i s 



(3.57) 
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IV. DISCUSSION 

In this Section, we examine more closely our results derived above, and compare them to the analysis in Ref. |45j . 

A. Free motor 

Let us first examine the simple motor properties in the absence of a cargo, i.e., for / = 0. Here the expression in 
Eq. I|3.41fl is at variance with the heuristic expression given by Mogilner |4q. 

9r + 9 s 2 

which is solely based on the reaction rates, neglecting the internal dynamics of the motor. To compare the expression 
for the present model, Eq. (|3.41[l we introduce the dimensionless parameters 



9r_ 

is'' 



(4.2) 



and 



q r = 9 -, (4.3) 

lr 

which express the ratio between the spring relaxation times, 77 1 and 7" 1 , and the residence times Eqs. <|3.20[l in 
states S and R, respectively. In terms of these parameters we obtain 

<«£> - («x)m x \ ■ (4.4) 

Cp + C» 1 + is + q r 

The correction factor to the heuristic velocity given by Mogilner et al. is clearly smaller than 1, but approaches 1 in 
the limit of £ p 3> Cu an d Qs <C 1, Qr *^ 1) which are exactly the conditions under which expression l|4.1|l was derived. 

We note that the velocity vanishes for ( p = Cv- In this case the attachment to the track has no effect on the friction 
and there is no motion. In the limit of large protein friction compared to the viscous drag coefficient, C, p 3> Cv, but 
q s ~ 1 and/or q r ~ 1 each conformational cycle does not yield a full step of length AL/2 due to incomplete spring 
relaxation, and the average velocity is reduced. In the limit of either q s 3> 1 or q r 3> 1 the motor comes to rest, as 
the relaxation rate or the hydrolysis rate become too large for the spring to change its average length. The motor 
would also function under conditions C > ( p or L s > L r , it would just move in the opposite direction. 

Inserting the characteristic biological numbers from Ref. [45| we have ( p /( v ~ 50, q r w 0.2, and q s w 2 • 10 -3 , and 
the correction factor takes a value of about 0.8. This corresponds to a average velocity of (v x }° ~ 4 x 10 3 nm/s. 
However, under different conditions, the discrepancy between the heuristic result (|4.1|) and the exact quantity 1)4. 4|) 
may become more significant. 

In Figs. and |3 we show the dependence of the load- free velocity on the hydrolysis rate g s and the protein friction 
( p (all other parameters fixed at the values of Ref. H^). Accordingly, with respect to these values the model motor 
velocity is close to optimum. The maximum in the (^-dependence shows the interplay between on and off-rates in 
the protein friction model, whereas the final plateau in the ^-dependence indicates the above-mentioned saturation, 
i.e., the motor still works for extremely large values of £ p , as long as q s and/or q r do not increase to high values, as 
well. We note that, as expected, the velocity goes to zero for vanishing hydrolysis rate, and when the protein friction 
approaches Cp — * Ct>- 

B. Motor carrying a load 

In the case of a load or cargo we proceed to discuss the expression for the motor mobility in Eq. (|3.57|) , which can 
be rewritten in the more convenient form 

= 2(g. r + g s )(P r k r + P s k s )(P r <; v + PsC P ) + k r k s (l + P r + PsCp/Cv) , 4 5) 
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In the further discussion of the mobility it is convenient to introduce the dimensionless parameters in Eqs. (|4.2|l and 
(|4.3|) . The expression l|4.5ll can then be reduced to the form 

(P r C„ + PsC P ){2q s + (1 + £)g r ) + (1 + P r )(v + P S (p 
M = 2C„(C»+Cp)(l + ? S +?r) ' (4 ' 6) 

Let us investigate this expression in some limiting cases: (i) In the absence of fluctuations, i.e., the case of a constant 
spring constant and rest length, in the relaxed state R, the protein friction £(t) is locked onto £ p , and we have P r = 1, 
P s = and q r = 0. By inspection of Eq. I|4.6[) we then obtain the mobility fi — 1/(C P + Cv), as discussed in Sec. 
Illb. (ii) Similarly, in the strained state S, the protein friction ((t) is locked onto ( v , P s = 1, P r = 0, and q s = 0, 
and we obtain the mobility \i = l/2£„. (iii) Finally, in the case ( p = ( v , we immediately find [i — l/2( v , as is also 
evident from the model equations (|2.1I2.2|) . Interpolating between the limiting cases (i) and (ii) above, we introduce 
the average mobility according to 

Mav = p r I h p s _L = + PrCv + PsCp _ ^ ^ 

Cp Cv ^Cv 2^(^p ~l~ Cu) 



and the mobility in Eq. I|4.6[) takes the form 



1 + ( 2g , + (1 + ^)q r ) 

At = Atav YT^TV r • (48) 

Inserting the characteristic biological numbers of Mogilner et al. l-V: Cp/Cu = 50, q r — 0.2, q s — 2 ■ 10~ 3 , and 
Pr ~ Ps ~ 0.5, we obtain the average mobility /i av « 2.6 x 10 5 nm/(s pN), while the correction factor in Eq. (|4.8|) is 
0.998, i.e., very close to 1. Hence, this gives rise to the ratio 



Mm 



13, (4.9) 



in comparison with the value estimated in Ref. |4E 



Mm = ~r- (4.10) 

4p 



The origin of the discrepancy between the present result and that of Ref. [45| is the slight difference between the 
models. In the strained state Ref. [45j operates with two characteristic times, that of S — > R conversion, i.e., the 
residence time Eq. (|3.20|) (t) s = l/g r , and the time of the restoration of bonds between the motor working head 
and the groove, which is much smaller. In the present model the two times are assumed equal, corresponding to the 
assumption that the spring relaxation S — > R is initiated when the working head becomes attached to the groove 
again. From a physical point of view this is equally possible, but implies that the motor spends much shorter times 
in state S than in state R, or g r 3> g s , i.e., P r ~ 1 and P s ~ 0, and the motor becomes much more volatile to the 
local force during these periods. 

From the mobility /i and the zero- load velocity (v®), we obtain the stall force 

(v°) 

/stall = ^Z. (4.11) 

M 

In Figs. 0| and we have on the same plots depicted the values of the mobility /i and the stall force / s t a ii versus the 
hydrolysis rate g s and the protein friction £„. We note that the stall force exhibits a maximum as function of g s close 
to the parameter values chosen in our calculations, whereas the mobility is close to saturation. Similarly, as function 
of £ p , the stall force is close to its maximum value, whereas the mobility does not change much within the chosen plot 
range (note that the ordinate does not reach the origin). In general, we observe that due to the particular dependence 
of /i on the model parameters, its value varies relatively weakly within large intervals for the individual parameter 
values. In Fig. Elwe have depicted the mobility and stall force as functions of the rate of relaxation g r . For large g r 
we obtain a stall force of the order pN, which is in the biological range. 



V. SUMMARY AND CONCLUSION 



In this paper we have by analytical means solved a molecular motor model proposed by Mogilner et al. |45|. This 
model represents a robotic motor solely based on an effective static friction interaction between motor and its support 
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(track). From the underlying Langevin equations, which represent the synchronized dichotomous processes of friction, 
effective spring constant and distance between the motor heads, we obtain explicit expressions for the load-free motor 
velocity, the mobility of the motor, and the stall force. Whereas the result for the load-free velocity produces a 
typical motor velocity of several (im/s for physiologically reasonable parameters, the exact solution overestimates the 
mobility, leading to a value for the stall force, which is roughly two orders of magnitude smaller than physiological 
values and significantly smaller than the simulations results reported in Ref . |45| . This variance is associated with a 
difference in the stochastic dynamics underlying the analysis and the dynamical processes implied in the numerical 
simulation. 

A likely explanation relies on the essential feature of the model, which is the decoupling of the dynamics of the 
motor protein-rail biopolymer interaction (chemical bonds forming and breaking) into a fast, detached process during 
energy consumption, and a slow relaxing process in one mechanochemical motor cycle. This purely stochastic picture 
leads to situations in which the motor detaches frequently, before its relaxing step is finished, and therefore the sub- 
cycles, which actually lend themselves to propulsion, are interrupted. Obviously, this leads to the underestimation of 
the stall force. In a real system, the fact that chemical bonds are established ensures that a full propulsion sub-cycle 
can be completed before dissociation takes place for the next loading of the internal motor 'spring' in parallel to 
hydrolysis. In comparison to the ratchet models in which the motor properties are represented by fluctuating between 
two different, periodic potentials, it appears that the latter rely on fewer parameters, and therefore their stall force 
can be adjusted better to actually observed values. 

We finally should like to emphasize that the obtained exact results allow for an exact and detailed study of the 
dependence of the motor characteristics on the various parameters without invoking numerical simulations. Additional 
features such as the low likelihood for detaching from the rail during the forward-motion, could be incorporated into 
the model and still be solved explicitly, using the solution schemes developed here. We therefore believe that this 
study leads to a better understanding of molecular motor models. 
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FIG. 1: Molecular motor model showing one mechanochemical cycle, during which internal fluctuations become directed 
('ratcheted') through protein friction. The conformational changes of the motor protein are represented by two states of an 
effective spring with rest lengths L s for the strained and L r for the relaxed states. 1 — > 2: Slow relaxation of the previously 
strained spring to assume the rest length L r \ due to protein friction, the white working head of the motor stays attached to 
the polar biopolymer (actin filament or microtubule), while the black idle head is free to move. 2-^3: As a consequence 
of ATP hydrolysis ('power stroke'), the spring contracts so quickly that the protein friction breaks down, and both heads 
symmetrically converge to assume the strained configuration with rest length L a . The distance covered per mechanochemical 
cycle is A = (L r - L s )/2. (Adopted from Ref. H^.) 
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FIG. 2: Dependence of the mean velocity (v®) on the hydrolysis rate g s , exhibiting a maximum at around 2.25 ■ 10 3 /s. The 
vertical dashed line shows the parameter value given in Ref. |4q|. For large values of g B , the velocity tends to zero. 




FIG. 3: Dependence of the mean velocity on the protein friction ( p , reaching a plateau for large values (f„ = 10 _6 pN 
s/nm). The protein friction £ p = 5 x 10~ 5 nm/(s pN) used in the calculations corresponds to the dimensionless value 50 in the 
plot, indicated by the dashed line. Note that = corresponds to Cp/Cv = !• 
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FIG. 4: Mobility \i and stall force / sta ii as a function of the hydrolysis rate g a . The vertical line marks the value from Ref. [45 
Note the different scales. 




FIG. 5: Mobility fi and stall force /stall as a function of the protein friction £ p (model value 5.0 ■ 10 5 nm/(s pN)). 
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FIG. 6: Mobility /i and stall force /stall as a function of the rate of relaxation g r . We observe that for large relaxation rate the 
stall force enters the biological range. 



